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4" Background Of T he Inventi on 



10J It has been recognized that there is a need to monitor the 

III 

P existence and effects of damage and corrosion in structural 

EQ 

materials. Aircraft components provide a much publicized 

Rl 

example, but the need exists in a variety of other structures, 
for example pipelines and building members. To detect defects, 

15 traditional modal analysis testing generally relies on 

structural excitation and vibration pickups. Structural 
excitation methods, in turn, usually employ either harmonic 
sweep or impulse excitation. The former method is more precise 
cpd can focus on resonant frequencies, while the latter is more 

20 expedient. Vibration pickups can measure displacement, velocity 
or acceleration. The use of self -conditioning accelerometers 
and laser velocimeters are known. Accelerometers allow 
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installation of sensor arrays that accurately and efficiently 
measure mode shapes, while lasers offer contactless measurements 
that are effective for low mass sensitive structures. 

In general, it is known to propagate electromagnetic 
signals through materials in order to measure anomalies in the 
materials. Ultrasonic transducers, for example, include 
piezoelectric crystals of a certain mass. The crystals are 
excited at a known resonant frequency so that the resonant 
crystal vibration generates an ultrasonic signal. Ultrasonic 
sensors, which are relatively bulky and expensive, are typically 
used in portable devices placed in communication with the 
structure of interest so that the crystals propagate ultrasonic 
waves into the material. Sensors detecting resulting signals in 
the material sense patterns that indicate anomalies in the 
materials such as cracks or corrosion. 

Summary Of The Invention 

It is an object of the present invention to provide an in- 
situ sensor system for structural material analysis. 

This and other objects may be achieved within a method of 
detecting a damage feature in a structure. 

A plurality of thin piezoelectric ceramic sensors are 
embedded on the structure. A first sensor is excited so that 
the first sensor produces a responsive signal in the structure. 
The responsive signal is received by a second sensor, and the 



presence or absence of at least one predetermined signal 
characteristic in the received responsive signal is determined. 
The predetermined signal characteristic is related to existence 
of the damage feature . 

5 The accompanying drawings, which are incorporated in and 

constitute a part of this specification, illustrate one or more 
embodiments of the invention and, together with the description, 

W serve to explain the principles of the invention. 

fen' 

SI 

jjjj Brief Description of the Drawings 

m 

1* A full and enabling disclosure of the present invention, 

fr 

a including the best mode thereof directed to one of ordinary 



fll skill in the art, is set forth in the present specification, 

n 

U) which makes reference to the appended drawings, in which: 

fas? 

R| Figure 1 is a cross-sectional schematic view of a thin 

15 piezoelectric sensor embedded on a structure for in-situ 

measurement in accordance with an embodiment of the present 
invention; 

Figure 2 is a cross-sectional diagrammatic illustration of 
a PZT wafer active sensor constrained by structural stiffness; 
20 Figure 3 is a graphical representation of a self -diagnostic 

performed in a sensor in accordance with an embodiment of the 
present invention; 
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Figure 4 is a graphical representation of an impedance 
spectrum measured by a sensor array in accordance with an 
embodiment of the present invention; 

Figures 5a and 5b are schematic illustrations of a sensor 
array in accordance with an embodiment of the present invention 

Figure 6 is an electrical schematic view of a sensor in 
accordance with an embodiment of the present invention; 

Figures 7a and 7b are schematic illustrations of a sensor 
array in accordance with an embodiment of the present invention 

Figures 8a-8d are graphical illustrations of elastic 
measurement waves generated by a transducer array in accordance 
with an embodiment of the present invention; 

Figures 9a- 9c are schematic illustrations of embedded 
sensors in accordance with an embodiment of the present 
invention; 

Figures 10a- lOd are graphical illustrations of measurement 
and response signals resulting from use of a transducer in 
accordance with an embodiment of the present invention; 

Figure 11 is a schematic illustration of a sensor array in 
accordance with an embodiment of the present invention in 
passive detection mode; and 

Figure 12 is a block diagram illustration of a sensor and 
monitoring system in accordance with one or more embodiments of 
the present invention. 



Repeat use of reference characters in the present 
specification and drawings is intended to represent same or 
analogous features or elements of the invention. 

Detailed Description of Preferred Embodiments 

Reference will now be made in detail to presently preferred 
embodiments of the invention, one or more examples of which are 
illustrated in the accompanying drawings. Each example is 
provided by way of explanation of the invention, not limitation 
of the invention. In fact, it will be apparent to those skilled 
in the art that modifications and variations can be made in the 
present invention without departing from the scope or spirit 
thereof. For instance, features illustrated or described as 
part of one embodiment may be used on another embodiment to 
yield a still further embodiment. Thus, it is intended that the 
present invention covers such modifications and variations as 
come within the scope of the appended claims and their 
equivalents. 

Generally, the present invention relates to embedding 
relatively small transducers on a structure in order to monitor 
the structure's structural health. Several examples are 
provided, namely, active electro-mechanical impedance sensing, 
active ultrasonic sensing, and passive sensing through 
interrogation of acoustic emissions and detection of elastic 
waves indicating impacts by a foreign body. A description of an 



electromechanical (E/M) impedance damage identification strategy 
is presented in the Appendix. Summarizing the identification 
strategy, the effect of a piezoelectric active sensor affixed to 
a structure is to apply a local strain parallel to the 
structure's surface that creates stationary elastic waves in the 
structure. The structure presents to the active sensor the 
drive-point impedance, Z J/r (co) = /(o/w e (co) + c e (Q))-/^(co)/co . Through the 
mechanical coupling between the PZT active sensor and the host 



ha? 

\l 

? £ structure, on one hand, and through the electro-mechanical 

ff-i 

l|£ transductance within the PZT active sensor, on the other hand, 

s the drive-point structural impedance is directly reflected into 

CSS, 

the effective electrical impedance as seen at the active sensor 
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C<? terminals, as shown in Figure 6. The apparent elect.ro- 

P 

rW mechanical impedance of the piezo-active sensor as coupled to 

15 the host structure is: 



/ 7 f,*\ 

Z(cq)= icoC 



1 K 3\ 



{G)) + Z str {G>) 

Where Z(o>) is the equivalent electro-mechanical admittance as 
seen at the PZT active sensor terminals, C is the zero-load 
capacitance of the PZT active sensors, K 3i is the electro- 
20 mechanical cross coupling coefficient of the PZT active sensor 

(k 31 =d u /yjs n s 33 ) , Z str is the impedance of the structure, and Z PZ t is 
the impedance of PZT active sensor. 
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The electro-mechanical impedance method is applied by- 
scanning a predetermined frequency range in the hundreds of kHz 
band and recording the complex impedance spectrum. By comparing 
the impedance spectra taken at various times during the service 
life of a structure, structural degradation and the appearance 
of incipient damage may be identified. The frequency should be 
high enough, however, that the signal wavelength is 
significantly smaller than the defect size. 

Referring to Figures 5a and 5b, consider an array of four 
active sensors. Each active sensor has its own sensing area 
characterized by a sensing radius and the corresponding sensing 
circles indicated at 18. Inside each sensing area, the sensor 
detection capability diminishes with the distance between the 
sensor and the damaged or corroded area. A damage feature that 
is placed in the sensor near field is expected to create a 
larger disturbance in the sensor response than a damage feature 
placed in the far field. Effective area coverage is insured 
when the sensing circles of several sensors overlap. The 
diagnostics of the adjacent structure is performed using the 
active (real) part of the E/M impedance. Incipient damage 
changes taking place in the structure are reflected in the 
drive-point structural impedance. The change in the structural 
drive-point impedance extensively affects the real part of the 



effective electro-mechanical impedance of the piezoelectric 
active sensor embedded on the structure. 

Figure 5a illustrates a structural crack placed in the 
sensing circle of the first active sensor. The presence of the 
crack modifies the structural field and effective drive-point 
structural impedance as seen by that sensor. At the same time, 
the crack also belongs to the sensing circle of the second 
sensor, but it is at the periphery of this circle. Thus, it is 
expected that the effective drive-point structural impedance as 
seen by the second sensor will also be affected but to a much 
lesser extent than that of the first sensor. The structural 
crack is outside the sensing circles of the third and fourth 
sensors, and their drive-point structural impedances will be 
almost unchanged. By virtue of the impedance equation above, 
changes in the drive-point structural impedance are directly 
reflected in the sensor's E/M impedance. Thus, the crack 
illustrated in Figure 5a is expected to strongly modify the E/M 
impedance of the first sensor, to slightly modify that of the 
second sensor and to leave the E/M impedances of the third and 
fourth sensors unchanged. 

Figure 5b illustrates a patch of corrosion damage placed in 
the sensing circle of the first active sensor. The corrosion 
damage also belongs, to a lesser extent, to the sensing circles 
of the second and fourth sensors. The effective drive-point 
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structural impedance seen by the first sensor will be strongly 
modified. The drive-point impedance seen by the second sensor 
will be modified to a lesser extent, and the impedance of the 
fourth sensor will be slightly modified. The drive-point 
impedance of the third sensor remains virtually unchanged. As 
indicated by the impedance equation above, these changes in 
drive-point structural impedance will be directly reflected in 
the sensors' E/M impedances. Thus, the corrosion damage 
strongly modifies the E/M impedance of the first sensor, 
somewhat modifies the E/M impedance of the second sensor, 
slightly modifies the E/M impedance of the fourth sensor and 
leaves the E/M impedance of the third sensor unchanged. 

The discussion above focuses on the use of E/M impedance 
sensing to detect defect or damage features in a structure. It 
is also possible to use active piezoelectric sensors to 
propagate elastic waves through the structure so that signal 
reflections and scattering effects may be used to detect damage 
features such as cracks and corroded areas. Although it should 
be understood that the present method may include the use of a 
single sensor, the embodiment described herein includes the use 
of a sensor array. In general, the sensors in a given array are 
clustered in a predetermined positional arrangement and are 
sequentially excited in round-robin fashion so that each sensor 
emits an ultrasonic acoustic wave into the structure. The 



remaining sensors in the array, as well as the originating 
sensor, monitor the emitted signals and, from the received 
responses, may determine the existence and position of a damage 
feature. Excitation may be at a fixed frequency, a frequency 
5 burst or a frequency sweep. The particular frequency band may 
be selected consistently with the size of the feature to be 
identified. 

n 

:E Figures 7a and 7b, illustrate an array of four active 

\J 

sensors in an ultrasonic sensing mode. Piezoelectric active 

m 

sensors can act as both sensors and actuators . At each 
5 measurement, one of the four active sensors is both an actuator 
ftj and a sensor, while the remaining three act only as sensors. 

q 

CO The actuator function moves round-robin among the four sensors. 
W Assuming the first sensor is the actuator, it generates elastic 
15 waves that propagate through the material and are sensed at the 
second, third and fourth active sensors. The properties of 
these waves are affected by the presence of impact damage or 
material defects and can be interpreted to yield damage location 
and amplitude. This method can be applied, for example, to 
20 detect damage arising from cracks and corrosion. 

Referring to Figure 7a, it is desirable to characterize 
crack damage in terms of its location (x c , y c ) , its size and its 
orientation (a,0) . When a crack is present in the wave path, a 
signal generated by the actuator piezoelectric sensor will be 
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deflected, reflected and transmitted at the crack. The relative 
proportions among deflection, reflection and transmission vary 
with damage, size and orientation. In Figure 7a, the second and 
fourth active sensors receive both direct and deflected wave 
signals from the actuator (i.e. the first active sensor). The 
first sensor also acts as a receptor and detects a wave 
reflected from the crack. The third active sensor receives a 
transmitted wave having an amplitude that is a function of the 
damage size. Thus, a matrix of information in terms of event 
arrival time and amplitude may be determined from the 
information detected by the four sensors when the first sensor 
is the actuator. Thereafter, a round-robin procedure is imposed 
in which the second, third, and fourth sensors in turn become 
the actuators. An information matrix is obtained at each such 
measurement. The damage, location, size and orientation may be 
derived using conventional linear algebra or neural network 
algorithms. Generally, the use of multiple round-robin 
measurements provides more information than needed. This over- 
specification, however, mitigates data error. 

Figure 7b illustrates a corrosion damage example. 
Corrosion damage is multi -dimensional and can cover a wide area. 
In the example shown in Figure 7b, the first active sensor 
generates elastic waves that propagate through the material and 
are sensed at the second, third, and fourth sensors. The 



damaged material affects the speed and attenuation of the waves. 
As with crack damage, it is desirable to characterize corrosion 
damage in terms of location, size and orientation. Orientation 
may be defined in any suitable manner, for example the major and 
minor axes of a damage ellipse and axis inclination. In the 
example, elastic waves received by the second and fourth sensors 
do not engage the corroded area and are, therefore, unaffected. 
Waves received by the third sensor travel through the damaged 
material and are therefore modified in speed and attenuation. 
The material degradation may be so advanced that no detectable 
waves are received at the third sensor. As each of the four 
sensors emits elastic waves in round-robin fashion, an 
information matrix is stored so that the location, size and 
orientation of the damaged area are determined. Again, the 
solution may be derived from conventional linear algebra 
solutions or neural networks, and the use of multiple 
measurements mitigates errors. 

The types of elastic waves used in the damage 
identification process may vary, for example among constant 
amplitude sine waves, bursts, sweeps and impulses. Four wave- 
type examples are provided in Figures 8a - 8d. The constant 
amplitude sine wave in Figure 8a is the simplest waveform. 
Because excitation frequency should be matched with the 
structural characteristics of the damage area in order to excite 
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a structural resonance, however, use of the constant amplitude 
wave may require multiple emissions at varying frequencies to 
test for all damage types of interest. Accordingly, a frequency 
sweep as in Figure 8b may be desirable in that it permits the 
5 excitation of more than one frequency in the same measurement. 
An impulse signal as in Figure 8c permits the excitation of a 
wide frequency spectrum. A frequency burst as in Figure 8d may 

i ' 

n 

j£j[. be desirable in that it includes a dominant frequency that can 

fll be tuned to particular structural and flaw- size requirements. 

yi 

15; Since its frequency content is known, a burst wave can be 

s readily detected and filtered from background noise, and its 

Q 

fll limited duration facilitates the identification and analysis of 

o 

ft? burst reflection from defects and boundaries. 

i» 3 3 

^ Experimental Results 

15 Lamb-wave propagation experiments were conducted on aging 

aircraft panels using a number of PZT active sensors affixed at 
various locations. Experiments were performed to verify the 
wave propagation properties and to _ identify reflections due to 
intrinsic construction features of the panels, for example, 

20 rivets, splice joints, etc. The illustrations in Figures 9a - 

9c provide a crack detection example. Referring to Figure 9a, a 
series of rivets 20 attach to panel sections 22 and 24 so that a 
side edge 23 of panel 24 faces toward the right in the view 
shown in Figure 9a. A piezoelectric active sensor 26 is 
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attached by adhesive to panel 24 at a distance of 200mm from 
rivets 20. A more complicated array of rivets is provided on 
panel 24 in Figure 9b. The panel shown in Figure 9c includes 
the same rivet arrangement, with the addition of a simulated 

5 crack (12.75mm EDM hairline slit) 28 starting from the right- 
most rivet at the top horizontal row. That is, Figure 9c 
illustrates the damaged panel. Experiments were also conducted 

p. on a panel (not shown) having no rivets or cracks. 

V. 5 

fjj In the clear panel (not shown) and the panels of Figures 9b 

m 

1$* and 9c, PZT sensors 26 emitted an approximately 25mV continuous 

sets 

s frequency wave. Figure 10a shows the emitted and received 

FIJ signals on the undamaged panel without rivets (not shown) . The 

W graph shows the emitted signal, centered at around 5.3 micro- 

fll 

seconds, and multiple reflections from the panel edges. The 
15 reflections start to arrive at approximately 60 micro-seconds. 

Figure 10b illustrates the signals for the panel shown in Figure 
9b. The signal features reflection from the rivets and multiple 
reflections from the panel edges. The reflection from the 
rivets arrives at approximately 42 micro- seconds , indicating an 
20 approximate time-of -flight having the value TOF = 37 micro- 
seconds. This TOF is consistent with traveling 200mm at a group 
velocity of approximately 5.4km/s. Figure 10c shows the signal 
recorded on the damaged panel of Figure 9c. The signal response 
is similar to that shown in Figure 10b, but is stronger at the 
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42 micro- second position. Subtracting the signal of Figure 10b 
from that of Figure 10c, as shown in Figure lOd, identifies the 
effect of the crack. That is, the strong wave pack centered on 
42 micro-seconds illustrates the response from the crack. 

Accordingly, where an array of sensors is disposed on or in 
a panel, and where signals are initially measured in the 
undamaged structure to obtain baseline measurements, the 
baseline measurements may be subtracted from subsequent 



measurements to identify the existence of damage features. The 

m 

LCr: time-of -flight determines the distance from the respective 



sensor. By measuring distances from multiple sensors in the 
array, the location of a damage feature with respect to the 
sensors may be determined. Moreover, the use of multiple 
measurements, sequentially activating the individual sensors in 
15 the array in round-robin fashion, provide for precise location. 

The examples provided above, i.e., impedancp^sensing and 
ultrasonic sensing, rely on active sensors^to measure structural 
characteristics. The present invention may also, however, be 
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^ised in passive modes to detect^structural damage. In certain 



20 v 0 f there embodiments, a plurality of sensors is disposed in a 
predetermined orientation relative to each other and at known 
positions on ther structure . The sensor outputs are monitored 
intermittently or continuously, even though the sensors may not 
be engaged in either of the active measurement procedures 
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described above. Damage events may 




identified through the 



reception of stress waves generated in the structure through 



impacts or other material disruptions. Certain waves may, for 
example, indicate an occurrence of a low- velocity impact. The 
sensor may also, hdwever, detect acoustic emission signals that 
indicate dam^e has occurred. By determining and recording the 
location/and time of damage events, a record may be compiled to 
preeact the structure's remaining operative life. 

Figure 11 illustrates a damage event taking place in a zone 
monitored by four piezoelectric transducers acting as passive 
sensors. In this example, the damage under consideration is a 
low-velocity impact. The impact location and intensity are . 
determined from elastic wave signals received by the four 
sensors. The arrival time at each sensor is proportional to the 
distance between the sensor and the impact location. 
Accordingly, a simple triangulation is used to determine the 
impact's position with respect to the sensors. This information 
is provided to a computer (not shown) that time-stamps the 
information so that a damage event record is kept. With four 
sensors, more data is available than necessary. Thus, an 
optimal solution algorithm, for example least squares, is used. 
Damage intensity is determined simultaneously by inclusion in 
the variable list. A neural network algorithm may be applied 
instead of explicit calculation. 




Information-rich data received from the sensors is 
preferably processed such that only relevant content is 
retained. The elimination of spurious noise and far-field 
disturbances, and the identification of useful information 
5 directly related to structural damage, is addressed through 
appropriate signal processing methodologies. Mathematically, 
the determination of the structure's physical condition using 

n 

sensor measurements is a nonlinear inverse problem. Signal 

Xtttf 

•~ a 

f|i processing and interpretation methods may be application- 
LO;; specific (model -based) or generic (non-model) . Conventional 

b (for example, Fourier analysis) and advanced (for example, 

□ 

W Wavelet analysis and digital filters) algorithms available in 

9 

jjjjj specialized software packages (for example, MATLAB, LABVIEW, 

ha? 

iy MATHCAD, or AUTOSIGNAL) may be used, as well as neural networks 

15 and expert systems. 

In operation, transducers in a measurement array execute 
self -diagnostics (as described in the Appendix) prior to damage 
detection measurement cycles. Temperature and humidity are the 
primary environmental changes that affect sensor calibration. 

20 Generally, these changes may affect the sensor itself and/or its 
adhesive interface with the structure. Although the values of 
the piezoelectric coefficients vary with temperature, the 
general aspect of the frequency response curve is only shifted. 
Thus, by calibrating the frequency shift over a variety of 
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temperatures, the system may compensate for temperature changes, 
thereby maintaining the general nature of the E/M impedance 
spectrum. The sensor/structure adhesion may be affected by both 
humidity and temperature, and the combination of humidity and 
temperature change may be especially detrimental to the adhesive 
interface. Thus, the self -diagnostic described in the Appendix 
with respect to Figure 3 may be employed to identify disbonded 
sensors . 

Each sensor may be formed as a small square or disk of thin 
piezoelectric ceramic material. In one embodiment, the length 
and width of each wafer is less than 13mm, +/-0 . 13mm, and 
between 0.2mm and 0.49mm thick, + /-0. 025mm, respectively. 
Electrodes are attached to the sensors, which communicate with a 
personal computer, other computing devices and/or other system 
components, through hard- wire or wireless connections. Sensors 
may be embedded on structures by, for example, insertion into a 
composite material during the layering process so that the 
sensor is actually within the structure or by adhesive 
attachment to the structure's outer surface as discussed above. 

Figure 12 provides an overview of an embodiment of the 
present invention. System 3 0 includes one or more structures 32 
having one or more transducer arrays disposed thereon as 
described above. For ultrasonic acoustic wave sensing, an 
ultra-high frequency high-energy tone-burst source 34 excites 
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individual sensors within each transducer array in round- robin 
fashion through a programmable wide-band bi-directional twenty- 
channel scanner/selector 36. These components communicate with 
a personal computer 38 through an ultra-fast data acquisition 
module 40. Computer 38 also acquires measurement data from the 
transducer array through data acquisition module 40. 

In wireless applications, the system may include a wireless 



W communication data-burst transmitter/receiver server 42 that 

„•« excites and reads the transducer array through a plurality of 

s II 

fs'S 

l¥ data-concentrators 44. That is, assuming there are a plurality 



3 of transducer arrays disposed about a structure, each transducer 

M 

. ffj is attached to an antenna that communicates with a data- 

p.% 

U) concentrator 44 dedicated to the transducer's array. The one or 
I'M more data-concentrators communicate with transmitter/receiver 
15 server 42 which, in turn, wirelessly communicates with an 

antenna associated with computer 38 or acquisition module 40. 

For impedance sensing, an adjustable-voltage-power gain- 
phase impedance analyzer 46 excites the t^nsducers, which 
output measurement information bapk'to analyzer 46. Analyzer 4 6 



2§)S includes software algarithp^rs to analyze the sensor data to 

determine the locatiprfs and orientation of damage features as 
described abaye\ Alternatively, analyzer 46 may forward the 
data to^omputer 38 for analysis. Generally, computer 38 houses 
tia^system's software components, the operation of which is 
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generally described above and which may include non- destructive 
evaluation and imaging software package 48, monitoring 



^^diagnostics and analysis software package 52 and/or artificial 

intelligence, neura/- network and data mining software 54. 

5 While one or more preferred embodiments of the invention 

have been described above, it should be understood that any and 
all equivalent realizations of the present invention are 

S included within the scope and spirit thereof. The embodiments 

VI 

flj depicted are presented by example only and are not intended as 

m 

l|g limitations upon the present invention. Thus, it should be 

BBS 

s understood by those of ordinary skill in this art that the 

0 

present invention is not limited to these embodiments since 
modifications can be made. Therefore, it is contemplated that 
any and all such embodiments are included in the present 
15 invention as may fall within the literal or equivalent scope of 
the appended claims. 



a 
m 
□ 
ru 



20 



it 
o 

rsi 
i 



E 



APPENDIX 

Because of their intrinsic electro-mechanical (E/M) 
coupling, commercially available low-cost piezoceramics can be 
used as both sensors and actuators . The frequency bandwidth of 
these devices is orders of magnitude larger than that of 
conventional shakers and impact hammers. Small piezoelectric 
ceramic wafers can be permanently attached to structural 
surfaces to form sensor and actuator arrays that permit 
effective modal identification in a wide frequency band. 
Assuming a single degree of freedom, the electrical admittance, 



□ as measured at the terminals of the piezoelectric (PZT) wafer 

sp 

ti attached to a structure, is as follows: 

5 

flj y(«)-tocfi-4 - z "<»> 1, (D 

I Z str (G>) + Z A (a>)) 

15 where C is the electrical capacitance of the PZT active sensor, 
Zstr (co) is the one-dimensional structural impedance as seen by 
the sensor, and Z A (co) is the quasi-static impedance of the 
sensor. At coupled- system resonance, the real part of the E/M 
admittance has a distinct peak. However, due to the additional 

20 stiffness contributed by the PZT wafer, the system's natural 

frequency may shift, although this effect diminishes greatly as 
the difference in size between the wafer and the structure with 
which it is used increases. 
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The E/M impedance/admittance response may be used to 
measure the structural dynamics of a material . Referring to 
Figure 1, consider a one -dimensional structure (for example a 
beam) 10 upon which is attached a ceramic PZT active sensor 12 

5 having a length l a/ at a position on beam 10 between x a and x a + 
l a . Upon activation, the PZT active sensor expands by zpzt. This 

ljl generates a reaction force F PZT from the beam onto sensor 12 and 

Q an equal and opposite force from the sensor onto the beam. This 

SI 

fll force excites the beam. At the neutral axis, the effect is felt 

m 

1^ as an axial force excitation, N PZ t, and a bending moment 

jL excitation, M PZT . As the active sensor is electrically excited 

i=J with a high-frequency harmonic signal, it induces elastic waves 

Hi into the beam structure. The elastic waves travel sideways into 

ru 

the beam structure, setting the beam into oscillation. In a 
15 steady-state regime, the structure oscillates at the PZT 

excitation frequency. The reaction force per unit displacement 
(dynamic stiffness) presented by the structure to the PZT 
depends on the internal state of the structure, on the 
excitation frequency, and on the boundary conditions: 

20 A«r(®)/fi«r(»)f (2) 

where u PZT (a>) is the displacement amplitude at frequency co, F PZT {u>) 

is the reaction force, and k str {®) is the dynamic stiffness. The 
symbol A signifies amplitude. Since the size of the PZT is very 
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small with respect to the size of the structure, equation (2) 
represents a point -wise structural stiffness. 

Consider a PZT wafer of length l a , thickness t a , and width 
b a , undergoing longitudinal expansion, u lf induced by the 
thickness polarization electric field, E 3 . The electric field is 
produced by the application of a harmonic voltage V(t) = Ke lCOt 
between the top and bottom surfaces (electrodes) . The resulting 
electric field, E = V/t, is assumed uniform with x x \-dEldx y =0 ) . 
The length, width, and thickness are assumed to have widely 
separated values ( t a «b a «l a ) such that the length, width, and 
thickness motions are practically uncoupled. 

Relevant equations regarding the PZT material are: 

S^sfc+d^ (3) 

D 3 = d 3] T ] +s 3 T 3 E 3 (4) 

where Si is the strain, T ± is the stress, D 3 is the electrical 
displacement (charge per unit area) , j/j is the mechanical 
compliance at zero field, e 3 r is the dielectric constant at zero 
20 stress, d 3i is the induced strain coefficient, i.e., mechanical 
strain per unit electric field. 

Since the piezoelectric active sensor is ultimately bonded 
to the structure, an elastically constrained sensor is 
considered. As a result of this assumption, the active sensor 

23 



is constrained by structural stiffness, k stri as shown in Fig. 2 
For symmetry, the structural stiffness has been split into two 
end components, each of size 2k str . The boundary conditions 
applied at both ends connect the resultant of internal stresses 
with the spring reaction force, i.e., 

t;(K)V.~2^,(K) (5) 

iH'.)v.=v,H'.) 

Substitution of Eq. (5) into Eq. (3) gives: 



(6) 



w ;(-i/ fl )^ fl =2^ r ^ Wl (4/ a ) + ^ 1 £ 3 

Introducing the quasi-static stiffness of the PZT active sensor, 



and the stiffness ratio 



r = Ac (8) 

k PZT 



Now Eq. (6) can be re -arranged in the form: 



/2 



Using Newton's law of motion, r/ = pu\ , and the strain- 
displacement relation, Eq. (3) yields the axial waves 
equation: 

",=c o y (10) 
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where { m y=d(ydt , and ()'=d()/ar, while c fl 2 = l/ps* is the sound speed. The 
general solution of Eq. (10) is: 

Wi(jc,0 = M^jOe"' where Wj(jc) = (Cj sinyx + C 2 cos/x) (11) 

The variable y = co/Co is the wave number, and O signifies .the 
harmonic motion amplitude. The constants C x and C 2 are to be 
determined from the boundary conditions. 

Substitution of the general solution (Eq. 11) yields the 
following linear system in Ci and C 2 : 



((p coscp + r sin cp)C, -(<psin(p -rcos(p)C 2 =\u !SA 
(<p cos <p + r sin (p)C } +(^sin^-rcos^)C 2 - \u lSA 



(12) 



fee* 

fy where u /SA =d 2] E 2 l a and <p=yY/ fl . Upon solution, 



> E 

W3B- 



IS" 



sinyjc 



2"/SA : 



?cos(p + rsin^ 



(13) 



Equation (4) can be re -expressed as: 



D 3 =-T(u'i-d i A)+eLE 3 = 4*3 



(14) 



where K,*, = </,*, is the E/M coupling factor. Integration of Eq. 
20 (14) yields the charge: 



2 2 ' |_ Vfl "31^3 ~2 / 



(15) 



where tl 2 b a tjt a = C is the conventional capacitance of the 



25 piezoelectric wafer. For harmonic motion, / = io> Q . Recalling the 
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expressions y = //^and Zsif 1 for the electric admittance and 
impedance, the admittance and impedance expressions for a PZT 
active sensor constrained by the structural substrate with an 
equivalent stiffness ratio r are: 



Y = iG>C 



1-K 2 
1 K 31 



1-- 



1 



<pcot<p 



4 



z = 



1 



\(oC 



i-k 3 2 ,|i- 



fi — i — T 

^ <pcot<p + r) 



(16) 



In this equation, the structural stiffness ratio, r, is 
additive to the sensor resonance term, (pcotq> . When the PZT 
active sensor is used in a frequency sweep, the apparent 
structural stiffness, k str , will vary with frequency, going 
through zero at structural resonances and through extreme values 
at structural anti -resonances . Eqs . (16) infer that both 
structural resonances and sensor resonances will be reflected in 
the admittance and impedance frequency spectra. 

Analysis of the asymptotic behavior of Eqs. (16) reveals 
important facts. , As r -» 0, i.e., vibrations of a free- free 



sensor, Y /ree = m-C 



1-K3.M- 



r. — i — t 

^ <pcot<p) 



As r becomes infinitely large, 



20 i.e., vibrations of a clamped sensor, Ydamped m C [ ] K «] . 0 n the 
other hand, as yla -» 0 (i.e., quasi-static sensor conditions) 

Y = ^. C \\^"], z = ^\\^"Y (17) 

L 1 + rJ \6)C\_ l + rj 
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The expressions contained in Eqs . (16) bridge the gap 
between high-frequency sensor- focused analysis and low- frequency 
structure- focused analysis. The present results cover the 
complete frequency spectrum and encompass both structure and 
sensor dynamics. 

The effect of structural and sensor damping can be easily 
introduced in Eqs. (16) by the use of complex notations: 



where 7 is the frequency- dependent complex stiffness ratio that 
reflects the structural point-wise dynamics and the sensor 
dissipation mechanisms. 

The response of the structural substrate to the PZT 
excitation is deduced from the general theory of beam 
vibrations. The excitation forces and moments acting upon the 
beam structure are derived from the PZT force, F pzr ^F PZT e^ , using 
the beam cross-section geometry: 



The space-wise distribution of excitation bending moment 
and axial force are expressed using the Heaviside function, H(x 
- x a ) , defined as H(x - x a ) = 0 for x < x a , and H{x - x a ) =1 for 





(18) 



M m -F m ±, N a =F PZT (19) 



N e (x y t)=N a [-H{x-x a )^H(x-x a -I a )}e^ (20) 

M e {xj) = -M a [-H(x-x a ) + H(x-x a -l a )\e i * t (21) 

Equations (2 0) and (21) correspond to axial and flexural 
vibrations, respectively. Axial vibrations modes are usually of 
5 much larger frequency than flexural vibration modes, but their 
vibration frequencies are commensurable with those of the PZT 
active sensors. Thus, both axial and flexural vibrations are 

□ 

□ considered. 

fljl The equation of motion for axial vibrations is: 

m 

l4* pA • u(x,t) - EA - u\x t t) = N'Jxj) (22) 

=|=r 

w Assume modal expansion 

. fSJ 

Q «M = ZCAW- e "' (23) 

CO 
15? 

F1J where X n (x) are orthonormal mode shapes, i.e., j x m x n dx = S mn , with 
5 mn = l for m = n, and 0 otherwise. C n are the modal amplitudes. 
Substitution of Eq. (20) into (22) yields: 

pA • u(xj)- EA • u'(x,t) = N a [S(x -x a )-^S(x-x a -l a )] e^ (24) 

20 where 5 is Dirac's function. Since the mode shapes satisfy the 
free-vibration differential equation EA-XZ+®* -pA-X n = 0 , 
multiplication by X n (x) and integration over the length of the 
beam yields : 

i N„ 



C -^M [ "^ ) + ^ +/fl)] (25) 
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Hence, 

u[xj) = Kfd^ (26) 

For Euler-Bernoulli beams, the equation of motion under 
moment excitation is: 
5 pA ■ w(xj) - EI • vf(x 9 t) - -M" e (x f t) (27) 

Substitution of Eq. (21) into (27) yields: 

o 

S| p/l • *(*,/) + £/ • w"(x,0 = M a [S'(x -x a ) + S'(x - x a - /„ )] ■ e*»' (28) 

ru 
cn 

1M where 5' is the first derivative of Dirac's function (6' = //'). 

Hi 

W Assume modal expansion h<jc,/)= Y C m X H {xy^ , where X n (x) are the 

ru 

Hj orthonormal bending mode shapes, and N lf N 2 are the mode numbers 

□ 

FM enduring the frequency band of interest. Since the mode shapes 
satisfy the free-vibration differential equation ElX'„'=<d];pAX n , 
15 multiplication by X n (x) and integration over the length of the 
beam yields : 

1 M n 



c > 



^lx„(x)[bXx-x a )^x^x a ^l a )]dx (29) 



" 2 2 a 
(O n - (O pA 

Integration by parts and substitution into the modal 
expansion expression yields: 

2 o c »=- r r L n^KW + ^^ + 0] (30) 

Hence, 
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v ( x ,,)=-tL f -f'M+r-i*. +0 . e *» (3i) 

pA£i a>-„-0> 2 
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JSMfl. 

u 
N 

ry 



fr 
is la a 



r«2 

■bar 
fir.* 

as a 

EM 



To obtain the dynamic structural stiffness, k str , presented 
by the structure to the PZT, the elongation between the two 
points, A and B, connected to the PZT ends is first calculated. 
Simple kinematics gives the horizontal displacement of a generic 
point P placed on the surface of the beam: 

u,(i)-u(x)-±W(x) , (32) 

where u and w are the axial and bending displacements of the 
neutral axis. Letting P be A and B, and taking the difference, 
yields: 

Using Eqs . (19), (26), and (31), the amplitude of Eq. (33) 
becomes 



15 



_ 1 PZT 



pA 



•(f)? 



(34) 



where differentiation between axial and flexural vibrations 
frequencies and mode shapes was achieved by the use of n Uf co„ h , 

U„(x) and n w , co„ , W n (x) , respectively. Dividing Eq. (34) by F PZT 

yields the structural frequency response function (FRF) to the 
20 Single Input Single Output (SISO) excitation applied by the PZT 
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active sensor. This situation is similar to conventional modal 



fcsr 

O 



testing with the proviso that the PZT wafers are unobtrusive and 
permanently attached to the structure. For convenience, the 
axial and flexural components of the structural FRF are 
expressed separately, i.e., 

#.(«)-— :Z" 2— T7 2-» (35) 



H>) .jJ»Y z IStflihSwr (3S) 



Modal damping, was introduced to provide practical 

veridicality to the model. The FRF 7 s are additive, and the total 
FRF is simply 

//(©) = //„(©) + //„(©) (37) 

The SISO FRF is the same as the dynamic structural 
compliance, as seen by the PZT wafer active sensor. . The dynamic 
structural stiffness is the inverse of the structural 
15 compliance, i.e., 



4r 
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□ 

09 
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ru 



U PZT 



(38) 



For free-free beams axial and flexural components [31] : 
^W = ^cos(y^), A^^lFl , y„.=^-, <»„.=y„c, c = j|, n u = 1,2,... (39) 

20 ^W=^[coshY %JC + cosy^jr-o^(sinhy^+siny ( , .»r)]<D J ^=y^fl , a = ^> A „w ml, ^lK( x ) dx ' (40) 
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Numerical values of /-y^and a„ w for w tt .<5 can be found in 



Blevins, R.D., (1979) "Formulas for Natural Frequency and Mode 



Shape," page 108, Litton Educational Publishing Inc.; for 5<n w , 



i(2« + i)ti 

= -± '— and a =1 . 

/ 2 



5 Experimental Results 

The analytical model described above was used to perform 

q several numerical simulations that directly predict the E/M 

sj impedance and admittance signature at an active sensor's 

rn 

. <if 

Ctt terminals during structural identification. Subsequently, 

s 

ItF experiments were performed to verify these predictions. The 



simulation conditions identically represent specimens consisting 
of small steel beams (E = 200GPa, p = 7750 kg/m 3 ) of various 

5 

f|j thicknesses and widths. All beams were 100mm long with various 



widths, bi = 8mm (narrow beams) and b 2 = 19.6mm (wide beams). 



15 The nominal thickness of the specimen was hi = 2.6mm. Double 



thickness specimens, h 2 = 5.2mm, were created by gluing two 



specimens back- to-back. Thus, four beam types were used: 



narrow-thin, narrow- thick, wide- thin, and wide- thick. The 



comparison of wide and narrow beams was aimed at identifying the 



20 width effects in the frequency spectrum, while the change from 



double to single thickness was aimed at simulating the effect of 



corrosion (for traditional structures) and 



disbonding/delamination of adhesively bonded and composite 
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structures. All specimens were instrumented with thin 7mm 

square PZT active sensors (1 = 7mm, b = 7mm, and t = 0.22mm) 

placed at x a = 40mm from one end. 

The numerical simulation was performed using a MATH CAD - 
5 coded simulation based on modal analysis theory, assuming a 

damping coefficient of 1% for steel beams. Numerically exact 
^ expressions for the axial and flexural frequencies and mode 
p shapes were used. The simulation was performed over a modal 

HI subspace that incorporates all modal frequencies in the 

yl 

1# frequency bandwidth of interest. The theoretical analysis 

-fS • 

indicates that these frequencies should be identical with the 



9 

Q 



jjj: basic beam resonances. The "calc" columns of Table 1 show the 

m 

HI first six predicted resonances for axial and flexural 

rii 

vibrations. During the experiments, a Hewlett-Packard 4194A 
15 Impedance Analyzer was used to record the E/M impedance real 
part spectrum in the 1 -3 0kHz range. To approximate the free- 
free boundary conditions, the beams were supported on common 
packing foam. The beam natural frequencies were identified from 
the E/M impedance spectrum. The results are given in the "EXP" 
20 columns of Table 1. When the beam thickness was doubled, the 
frequencies also doubled. This is consistent with theoretical 
prediction. The error between theory and experiment, however, 
seems larger for the double thickness beam. 



33 



The experimental results indicate frequency clusters, which 
move to higher frequencies as beam width is reduced. The 
frequency clusters appear to be associated with width 
vibrations. Width vibrations are also influenced by thickness, 
5 i.e., they shift to higher frequencies as thickness increases. 
This is also noticeable in Table 1, which shows that the lowest 
cluster appeared for the single- thickness wide specimen and that 

P 

p the highest cluster appeared for double -thickness narrow beam. 

M 

ill For purposes of comparison, dynamic identification was also 

1$ attempted with conventional modal analysis methods. A small 

B5CI 

^ steel beam with dimensions identical to the "narrow thin'' 

Q 

■il? specimen in Table 1 was instrumented with two CEA-13-240UZ-120 

j= 

Jlj strain gauges connected in half bridge configuration to a P-3500 

ft! 

' !r strain indicator available from Measurements Group, Inc. The 
15 specimen was suspended in a free-free configuration and excited 
with a sharp impact. The resulting signal was collected with an 
HP 54601B digital oscilloscope and numerically processed on a 
personal computer. Standard signal analysis algorithms (FFT) 
were used to extract the frequency spectrum. The first natural 
20 frequency (1.387kHz) was clearly displayed. The second natural 
frequency (3.789kHz) could also be identified but with a much 
weaker amplitude. These results were consistent with 
theoretical values and the experimental results presented in 
Table 1, but the impact excitation method was not able to excite 
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other higher frequencies depicted in Table 1, most probably due 
to bandwidth limitations. 

In order to obtain consistent experimental results, a set 
of sensors was subject to a calibration procedure consisting of 
5 geometrical and electrical measurements and measurements of the 
intrinsic E/M impedance and admittance spectra of the PZT active 
sensors. The results show small dimension variations from 

jUjjj sensor to sensor. The mean and standard deviation values for 

Sj 

ry length/width and thickness were 6.95mm, +/- 0.5%, and 0.224mm, 



ms +/- 1.4%, respectively. Electrical capacitance was measured on 

HP 

5 an in-process quality check to be applied during each step of 

q 

HI sensor development and during the sensor installation process. 

O 

W Mean and standard deviation values were 3.276mF, +/- 3.8%. 

a 

Hi 

1 The intrinsic E/M impedance and admittance spectra of the 

15 PZT active sensors, before being attached to the structure, were 
measured with an HP 4194A Impedance Phase-Gain Analyzer. The 
PZT wafer was centered on a bold head and held in place with a 
probe tip. Thus, center clamping conditions were simulated, and 
the wafer could perform free vibrations while being tested. 
20 Mean values of 251kHz and standard deviations of +/- .2% were 
obtained. 
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Table 1 



Theoretical and experimental results for wide and 
narrow beams with single and double thickness 



# 


Beam #1 (narrow thin) 


Beam #2 (Narrow thick) 


Beam #3 (wide thin) 


Beam #4 (wide thick) 


Calc. 
KHz 


Exp. 
kHz 




Calc. 
kHz 


Exp. 
kHz 


A% 


Calc. 
kHz 


Exp. 
KHz 


A% 


Calc. 
kHz 


Exp. 
kHz 


A% 


1 


1.396 


1.390 


-0.4 


2.847 


2.812 


-1.2 


1 .390 


1.363 


-1.9 


2.790 


2 .777 


-0.5 


2 


3.850 


3.795 


-1.4 


7.847 


7.453 


-5.2 


3 .831 


3 .755 


-2 


7.689 


7.435 


-3.4 


3 


7.547 


7.4025 


-2 


15.383 


13.905 


-10.6 


7.510 


7.380 


-1.7 


15.074 


13.925 


-8.2 


4 


12.475 


12.140 


-2.7 




29.650 




12.414 


12.093 


-2.6 




21.825 




5 


18.635 


17.980 


-3.6 


25.430 


21.787 


-16.7 


18 .545 


17.965 


-3.2 


24.918 


22 .163 


-12.4 


6 




24 . 840 












24 .852 










7 


26.035 


26 .317 


1 


26.035 


26.157 


0.5 


26.022 


26.085 


0.2 


25.944 


26.100 


0.6 






Cluster 
175 kHz 






Cluster 
210 kHz 






Cluster 
35 kHz 






Cluster 
60 kHz 





The active sensors used in the experiments were very small 



and therefore did not significantly disturb the dynamic 
properties of the structure under consideration. Table 2 
presents the mass and stiffness for the sensor and structure. 
For purposes of comparison, the table also presents the mass of 
an accelerometer . 



Table 2 



Numerical illustration of the non- invasive 
properties of the piezoelectric wafer active sensors 



Sensor type 


Mass, g 


% of 
structural 
mass 


Stiffness, 
MN/m 


% Of 
structural 
stiffness 


Active sensor 


0.082 


0.5% 


15 


1.5% 


Structure 


16.4 


N/A 


1000 


N/A 


Accelerometer 
:352A10, PCB 
Piezot ronics 


0.7 


4.3% 


N/A 


N/A 



Table 2 illustrates the non-invasive properties of 



piezoelectric wafer active sensors. The mass and stiffness 
editions of the sensors are within about 1% (0.5% for mass, and 
1.5% for stiffness). 
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* 



As described in more detail above, the present invention 
contemplates embedding piezoelectric wafer transducers on 
mechanical structures for potentially long periods of time. 
Thus, it is desirable to intermittently monitor the integrity of 
5 a transducer or transducer array. Accordingly, a self- 
diagnostic may be employed in which the reactive part of the 
transducer's impedance is measured over a relatively broad 

jack 

n 

m frequency range. Piezoelectric active sensors are predominately 

hat* 

H 

rfj capacitive devices dominated by reactive impedance. Baseline 

[?? 

LG; impedance signatures can be used to identify defective active 

ic 

s !«r 

b sensors. Referring to Figure 3, for example, a reactive 

5 

»M . impedance spectrum 14 of a well-bonded PZT sensor is compared > 

9 

with an impedance spectrum 16 of a disbonded (free) sensor. The 
^ appearance of sensor free -vibration resonance and the 

15 disappearance of structural resonances provide un-ambiguous 

features that indicate the degree of the transducer's integrity. 

Piezoelectric active sensors, and the associated structural 
dynamic identification method described above based on electro- 
mechanical impedance response, is well-suited for small rigid 

2 0 machinery parts having natural frequencies in the kHz range. 
For example, Figure 4 illustrates the natural frequencies 
identified by two PZT active sensors disposed on an aircraft 
turbo-engine blade, one on the blade's main portion and one at 
the blade root. Following identification of such impedance 
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spectra from the sensors, the sensors may thereafter be 
periodically monitored for deviations in the spectra that 
indicate structural defects in the blade. 
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